Instability Crossover of Helical Shear-flow in Segregated Bose-Einstein Condensates 
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We theoretically study the instability of helical shear flows, in which one fluid component flows 
along the vortex core of the other, in phase-separated two-component Bose-Einstein condensates 
at zero temperature. The helical shear flows are hydrodynamically classifled into two regimes: (1) 
a helical vortex sheet, where the vorticity is localized on the cylindrical interface and the stability 
is described by an effective theory for ripple modes, and (2) a core-flow vortex with the vorticity 
distributed in the vicinity of the vortex core, where the instability phenomena are dominated only 
by the vortex-characteristic modes: Kelvin and varicose modes. The helical shear-flow instability 
shows remarkable competition among different types of instabilities in the crossover regime between 
the two regimes. 

PACS numbers: 67.85.Fg, 03.75.Kk, 47.37.-fq, 67.85.De 



I. INTRODUCTION 

Kelvin- Helmholtz instability (KHI), one of the most 
fundamental instabilities in fluid dynamics, occurs in the 
presence of shear flow between two immiscible fluids [l| . 
When the relative velocity between the two fluids is suf- 
ficiently large, a vortex sheet existing along the inter- 
face between the two becomes unstable and develops into 
characteristic roll- up patterns of eddies. Because of the 
universal applicability of fluid dynamics, KHI may occur 
not only in classical fluids but also in quantum fluids such 
as superfluids [3,l3|- However, the KHI in superfluids — 
namely, quantum KHI — can be distinct from that in a 
classical fluid since the dynamics is governed by macro- 
scopic quantum effects, superfluidity, and vortex quanti- 
zation. Recently, quantum KHI has been proposed in 
phase-separated two-component Bose-Einstein conden- 
sates (BECs) consisting of two distinguishable Bose par- 
ticles [3, Q. Quantum KHI is realized in the presence 
of a relative superfluid velocity between the two compo- 
nents. When the relative velocity exceeds a critical value, 
a flat interface between the two superfluids changes into 
characteristic sawtooth waves without energy dissipation, 
leading to the formation of single-quantized vortices from 
the peaks and troughs of the waves, which is quite differ- 
ent from the nonlinear development in KHI in classical 
fluids. 

Although quantum KHI is a quantum counterpart of 
KHI in classical fluids, exotic hydrodynamic instabilities 
without classical counterparts also offer interesting pos- 
sibilities Q. Such instabilities can occur as a direct re- 
sult of the quantum effects in superfluids: superfluidity 
and vortex quantization. For example, superfluidity — 
the disappearance of friction — makes it possible to real- 



ize a counterflow of two miscible fluid components, called 
countersuperflow, in two-component BECs. The insta- 
bility of countersuperflow leads to exotic instability phe- 
nomena [y, lZI and has been realized very recently in two- 
component BECs [g . Another possibility can be realized 
in the presence of quantized vortices. Kelvin wave (or 
Donnelly-Glaberson) instability [9|, |lO[, which causes a 
deformation of a vortex from a straight line into a helix, 
is one such fundamental type of instability discussed in 
different superfluids |ll| - [l3j . 
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FIG. 1: Schematic diagram of helical shear flow in phase- 
separated two-component BECs. The first component forms 
a quantized vortex and the second component flows along the 
vortex core. 

In this work, we theoretically study the instability 
of helical shear flow in phase-separated two-component 
BECs, in which one component flows along the core of a 
quantized vortex of the other component forming a cylin- 
drical interface, as is illustrated in Fig. [TJ This state is 
called helical shear flow since the trajectory along the 
local relative velocity between the two fluids is helical. 
Helical shear flow is a flow state stabilized by the invis- 
cid flows and vortex quantization in superfluids. Helical 
shear flows may naturally arise in turbulent or nonequi- 
librium states in segregated condensates with a large pop- 
ulation imbalance between the two components (e.g., see 
Refs. [ij-[l3)- ^ helical shear flow becomes dynami- 



cally unstable when the relative velocity along the vortex 
core exceeds a critical value. The instability phenomena 
change drastically depending on the rotational velocity 
and the linear velocity of the first and second compo- 
nents, respectively. In this paper, we develop the phase 
diagram and investigate the nonlinear development of the 
helical shear-flow instability. 

The paper is organized as follows. In Sec. II, we char- 
acterize the stationary helical shear flow from its vorticity 
distribution in the Gross-Pitaevskii (GP) model at zero 
temperature. Sections III and IV discuss the linear sta- 
bility of the helical shear flows and show an instability 
phase diagram on the basis of an effective theory and nu- 
merical analyses of the Bogoiubov theory. In Sec. IV, 
we demonstrate the nonlinear development of the heli- 
cal shear-flow instability by numerically solving the GP 
equation. Section V is devoted to the conclusions and 
discussion. 



II. STATIONARY HELICAL SHEAR FLOW 

Let us consider two-component BECs in a uniform sys- 
tem in the GP model at zero temperature [22]. Two- 
component BECs are well described by the complex order 
parameter ^^ (r , t) of the jth component. The mean-field 
Lagrangian is given by 



C= I dPx [Pi + Vi 



.gi2|*i|'|*2|'] 



(1) 



with 



where rnj and /ij arc the particle mass and the chemical 
potential of the jth component, respectively. The inter- 
action parameters g^-fe arc defined as gj^ = iirh^ajkirrij + 
mk)/mjmk (j, fc = 1, 2) with the s-wave scattering length 
Ujk between the jth and kth components. The parame- 
ters are assumed to satisfy the condition for phase sepa- 
ration, 512/^511522 > 1- In the following, we shall con- 
sider experimentally feasible situations, m = nii = 7712, 
9 = 911 = 522, and 512/5 = 1-2, e.g., weakly seg- 
regated condensates in Ref. [l^. Considering station- 
ary states with rotational symmetry along the z axis, 
*j(a;,t) = <i>j{x) = ^nj(r)e'('3j^+^^'^) with the cylin- 
drical coordinates x = {r,(f),z), one obtains the time- 
independent GP equations 



Yl 9jknk 

k 



Hj 



$, = 0. 



(3) 



Here, we used hj = -~^ {df. + r-^dr - Q] - L^r^). 
Because we are interested in the vortex state illustrated 
in Fig. [U we set without loss of generality 



where the first component has a linear velocity — paral- 
lel to the z axis and the second component has a ro- 
tational velocity — around the axis. In the follow- 
ing, we use the dimensionless parameters 7 = 512/5 

and ly = /^2/Mi with /i2 = 1^2 — n '^ ^^'^ ^^'^ scales 
^ = y^h?/{m^i), T = h/ni, and c = y^fii/m for length, 
time, and velocity, respectively. 
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FIG. 2: Radial profiles of the density of two components, 
|$jp, and vorticity iv in the stationary state with 7 = 0.95 
and V — 1.2. The radius ro of the interface, defined by the 
plane |$i| = |$2|, increases with \L\. The horizontal axis is 
the radial coordinate. The vertical axises on the left and right 
hand sides shows the density and the vorticity respectively. 

Figure [5] shows radial profiles of characteristic heli- 
cal shear flows in the stationary states with 1/ = 0.95. 
When the first component contains multiquantized vor- 
tices [L = 4 in Fig. [2], the centrifugal force resulting 
from the rotational flow in the first component enlarges 
the radius rg of the circular interface plane, defined by 
the plane |^i| = |^2|. The interface radius ro becomes 
smaller as L decreases. The interface is ill-defined for 
i = 1 in Fig. [2l^d), where the radius tq becomes compa- 
rable to the interface thickness ~ ^. 

To classify the stationary shear-flow states and char- 
acterize their hydrodynamic aspects, it is convenient to 
introduce the mass-current velocity v and the vorticity 
w, defined as 



1^^. [**(VvI/,)-(VvI/*)vI/,] 



E,"^jl*.f 



w 



V X w. 



(5) 
(6) 



(Qi, Q2, ii, i2) = (o, g, L, 0), 



(4) 



For L = 4 in Fig. [21 the velocity v changes sharply and 
the vorticity w = (0, w^, Wz) is localized around r = tq. 



where the vortex sheet forms a circular pipe from the 
hydrodynamic viewpoint. Since the trajectory along the 
direction of w is helical, this state may be called a helical 
vortex sheet. When the interface radius tq becomes close 
to the z axis (r = 0) for L = 1 in Fig. [2l the density 
71,2 (r = 0) at the center is suppressed so that the vorticity 
is distributed around r — 0, forming a linear vortex with 
a superflow along the vortex core. We call such a vortex 
state a core-flow vortex. 

The interface radius tq depends on the particle pop- 
ulations of the two components. Because the popula- 
tions of the first and second components are decreasing 
and increasing functions of v with L fixed, respectively, 
the radius rg is an increasing function of v. There ex- 
ists a critical ratio (t'min) of u, below which the popu- 
lation of the second component vanishes in the vortex 
core. In the limit of $2 ^^ for h' — > i^min, by ne- 
glecting the n2 terms in Eq. ^, $2 corresponds to the 
single-particle wave function of the ground state in an 
external potential gi2'^i- Consequently, the value ^'minMi 
equals the ground-state energy; e.g., for L = 1, we have 
i^min = aiV7 - 0.899 with m = faJir/O^ {r ^ 0). 
However, we found that the interface radius ro is larger 
than the system size for :/ — > 1 in the numerical simula- 
tions. This fact may be understood from the hydrody- 
namic viewpoint discussed in the next section: the small 
difference between the hydrostatic pressures on the two 
sides of the interface balance for sufficiently large tq . Al- 
though the instability discussed below also depends on 
the ratio v, we consider the instability for v = 0.95 in 
the following; this case is sufficient to show the charac- 
teristic behavior of the shear-flow instability. 



III. EFFECTIVE THEORY OF HELICAL KHI 

First; we perform a linear stability analysis of the he- 
lical vortex sheet. The effective theory of quantum KHI 
of a flat vortex sheet is applied to our case of a helical 
vortex sheet. In the effective theory, we consider only 
the radial shift 77 of the interface position and neglect the 
thickness. The position of the interface is represented by 
a curve r = rg + r]{(j), z, t). 

In this approximation, the Lagrangian [Eq. ([1])] is re- 
duced with the interface-tension coefficient a and the in- 
terface area S to 



/•CXD 

^Rff = / dz 

-00 

-aS. 



rdrVi + / rdrV2 

ra+ri Jo 



(7) 



The area S equals Sq = 2tt J dzr^ in the stationary sate 
with 77 = 0. For a small deformation, we have 



S^So 



dzd(j> 



Zro I 



(8) 



The variation of Eq. ([7]) with respect to t] gives 
-(ro + T])Vi{ro +1]) + (ro + 7])'P2{ro + ??) 



+ ai-l + ro[ -^d^-q + a/?7 n. = 



ro^ 



(9) 



This equation shows the pressure equilibrium on the in- 
terface in the stationary state. 



P2-P1 



a 



(10) 



where Pj is the value of Eq. ([2]) along the interface in the 
stationary state, called the hydrostatic pressure. The 
left- and right-hand sides of Eq. (|T0)) refer to the pres- 
sure difference across the interface and the tension result- 
ing from the curvature of the interface, respectively. In 
the Thomas-Fermi approximation neglecting the quan- 
tum pressure term oc (^rij)'^ /rij in Eq. ([2]), one obtains 



P. 



^gn'^Aro) with 77 1 



Ml 






/g and 772 ^ fJ-^g- 



The equilibrium radius ro is calculated from Eq. (|T0|) . 
For large ro, one can take the coefficient a in Eq. (fTU)) 
to be that for a flat interface, calculated by assuming 

Pi^P2 MM, 



aov^ 



(11) 



with the tension coefficient 

the external pressure y.\/2g. 
o (i^C'^/ro) , we obtain 



ao = "^Vt^I under 
Neglecting the term of 






^"(7-1) 



L2(i_ 



■^"(7-1) 



1-7/2 



(12) 



This equation holds for l/\/2 < ro/L^ < l/\/2(l - v) 
from the conditions 77i(ro) > and P2 > Pi- The ex- 
pression p2p is in good agreement with the numerical 
results of the interface radius of the helical vortex sheets 
in Fig. [5] [see also the insert in Fig. [3], and thus, we can 
compare the stability analysis demonstrated below with 
the numerical analysis presented in the next section. 

Linear stability is investigated by considering a small 
fluctuation, 



|vl/,(a;,t)| 



with 



[^ = 77j(r) +6nj{x,t), 


(13) 


t) ^ QjZ + Lj<f> + Sej{x,t), 


(14) 


f] sin(/0 + qz — ^i,qt), 


(15) 


Snj(r)sm{l(j) + qz — (jJi.qt), 


(16) 


S9j{r)cos{l(j) -j- qz -- '^i,qt), 


(17) 



50, 



where I and q are real. The fluctuation represents the 
ripple mode, which causes a ripple wave on the interface. 
The kinematic boundary condition is employed along the 
interface. 



dtV' 



Vt], 



(18) 



with the radial component [vj], of the local superfluid 
velocity V, = ^-^ [**(VvI/,) - {V^*)^,] of the jth 

component. The radial shift of the interface causes a 
density fluctuation along the interface, Sni = —rjdrrii = 
—rj i and 5n2 — —ridrn2 = 0. The boundary con- 

dition (|18p. together with the equation of continuity, 
dtSuj + V (rijdvj) + V {SrijVj) = 0, yields [V • 5vj]^^^^^ = 

with 5vj = '^^^^j- Since the fluctuation is consid- 
ered to be localized around the interface, we choose a 
form of the phase fluctuation as 



S0i{r) 



CiKi{\q\r), 
C2li{\q\r), 



(19) 
(20) 



where ci and C2 are constants and Ii and Ki are modified 
Bessel functions of the second and first kind, respectively. 
By linearizing Eq. (O together with Eqs. dTB]), (IT9l) . 
and (|^ . we obtain the dispersion relation 



J- /Kn ^ -i ' ii lyy-t V 



Wj, 



^1 -1-^2 



(21) 



with 



A = (Si+S2)|g| 



/ + a 



-S1S2 [q ■ VeY 



(22) 



Here, we used .=,1 



2mn2(ro)Ii(\q\ro) 



r,/ 



2mni(ro)_fi'i(|g|ro) 
Ki-i(\q\ra)+Ki + i{\q\ro 



-ni(ro)- 



(0,V»"o,'7), 



Il-l(\q\ro)+Ii + i(\q\ra) 

and Vr^ = m e^' L/rQ, —Q). The frequency a;;_q is real for 
A > but becomes complex for ^ < 0. For ^ < 0, 
the helical vortex sheet is dynamically unstable against 
some perturbation so that ripple modes are amplified ex- 
ponentially with time. For A > 0, the fiLOi^q may corre- 
spond to the energy quanta of ripple modes or ripplons. 
In the presence of energy dissipation, the helical vortex 
sheet is thcrmodynamically unstable when a ripple mode 
has negative energy with cj; ^ < 0. Here, we restrict the 
discussion to dynamic instability by neglecting thermo- 
dynamic instability, which is experimentally reasonable 
in atomic BEC systems. 

The dispersion relation ([?!]) has the generalized form 
of that of quantum KHI [j| and capillary instability [23] 
in phase-separated two-component BECs. For L = Q = 
I = 0, where there is no shear flow between the two com- 
ponents, the dispersion relation (PT|) reduces to those of 
the Plateau- Raylcigh (or capillary) instability [2j]. In 
general, a capillary flow, namely, a helical vortex sheet 
with L = 0, is dynamically unstable because A < with 
/ — — a/r-Q in Eq. (22), whereas a helical vortex sheet 
with L 7^ can be stabilized because of the centrifugal 
force in the presence of a vortex at r = 0. However, in the 
limit of To — )■ 00, the dispersion relation (PT|) reduces to 
that of the quantum KHI for a flat interface without ex- 
ternal potentials in Ref. |3|. In this sense, the instability 



for the case with a sufficiently large tq may be called the 
helical KHI. In the same manner as in those instabilities, 
our effective theory is based on the assumption that the 
fluctuations are well localized around the interface and 
the wavelength of the ripple modes are smaller than the 
thickness of the interface. 

The phase diagram for the linear stability is calculated 
from the dispersion relation (pij) . The frequency (PT|) of 
a ripple mode with / becomes complex when Q exceeds 
the critical value 



Qi = mm{Qi,q), 
q 



(23) 



with 



ii.q 



q\ 



,■ ^1 i-'^2 
^l'^2 



J + aiP/r^ + q^)] 



L\l\ 



We plotted the L dependence of Qi with 1/ ~ 0.95 and 
7 = 1.2 for < ; < 6 in Fig. H The helical shear flow 
becomes dynamically unstable when the velocity — Q of 
the second component exceeds the critical value 



Vc = —mm{Qi). 
ml 



(24) 



The critical velocity Vc decreases monotonically with L 
as seen in Fig. [3] The velocity Vc approaches zero asymp- 
totically for L — > 00 since then the interface radius ro in- 
creases with L as ~ £,L/y/l — v and then the stabilizing 
force / is asymptotic to zero. 
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FIG. 3: L dependence of Qi in Eq. ((231 for 1 < L < 4 with 
V = 0.95 and 7 = 1.2. The insert shows the L dependence of 
the interface radius ro in Eq. (|12p . 



IV. CROSSOVER FROM HELICAL KHI TO 
CORE-FLOW INSTABILITY 

When the interface radius rg becomes small and the 
helical vortex sheet turns into a core-flow vortex, the rip- 
ple modes are ill-deflned, and thus, the effective theory 



can no longer be applied. To determine the stability in 
the crossover regime between the helical vortex sheets 
(ro ^ oo) and the core-flow vortices (ro — > 0), wc demon- 
strate numerically the regime's linear stability based on 
Bogoliubov theory [221 • 

By introducing a perturbation of the order parameters, 

S^j := g*(L.0+Q.^) (j^^ -~V*), (25) 

and linearizing the time-dependent GP equations, 



ihdt^ 



't^j 



2m 



■V' + E5j^l*jl'-^0*^' (26) 



around the stationary states ^j = $j with Uj = 
Mj(r)e^('J^+''^-"*) and Vj = „j.(r)e*(9^+''^-'^*), we obtain 
the Bogoliubov-de Gennes (BdG) equations, 

huju = a^Hu, (27) 

where we used u = {ui,U2,vi,V2)'^ , 13 = 
diag(l, 1,-1,-1), 



fh^ + Gil 0x2 



Gil 



n 



Q12 

Gil 
V G12 



G12 \ 
G22 



G12 h^ + Gil G12 



(28) 



G22/ 



with Gjk = gjky^njUk, and 

1 il±L,r 



2m 



'J 



d^ + -dr~'-^^^~iq±Q,r 



-i-h 



Y,93knk. 

k 



(29) 



An eigenvalue uj of Eq. (P7| represents the frequency of an 
elementary excitation in the helical shear-flow state. Be- 
cause a solution (w, (7, /, Uj,Vj) has its conjugate solution 
(— w*, — g, — Z, V*, u*), we present here, only the results for 
q> Q without loss of generality. 

Figure |3] shows the distribution of the imaginary part 
Qoj > for several values of I obtained by numerically 
solving the BdG equations (p7)) . The effective theory well 
explains the imaginary part of the the numerical results 
for all I in the limit ro — > 00 (see, e.g., the plot for L = 10 
in Fig. 2]). However, in the crossover regime where the 
interface radius tq becomes comparable to the interface 
thickness ~ S,, the imaginary part for the specific mode 
|/| = L > 1 of the numerical results deviates from that 
of the analytic result compared to the other modes, as is 
seen in Fig. |4]for L > 1. The deviation is related to a 
specific behavior in the spatial distribution of the |Z| = L 
mode. We can see that the distribution |(5^ip has a 
finite value at r = only for |^| = L (see Fig. \5^. This 
is because the amplitude of ui (vi) can survive at r = 
for I = — L {I = L) with the third term in Eq. (P^]) being 
zero. However, the modes with \l\ j^ L are distributed 
locally around the interface r = rg in Fig. [51 and then the 
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FIG. 4: The imaginary part 9a; > for several values of 
I obtained by the analytical (curves) and numerical (points) 
calculations based on the effective theory and the Bogoliubov 
theory, respectively. We set the parameter Q^ = 1.23 for 
L = 1 and Q^ = 0.74 for L = 2, 3, 4 with v = 0.95 and 
7 = 1.2. 
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FIG. 5: The radial distribution of \5'^i{r,(j>,z,t)\'^ of the 
first component for <j) — z — t = obtained by numerically 
solving the BdG equations (|27|l . The splitting modes (1Z| = L) 
have a finite value at r = while the modes with |Z| 7^ L are 
distributed locally around the interface r — ro- We set the 
parameter Q^ = 1.23 for L = 1 and Q£, = 0.74 for L = 2, 3, 4 
with u = 0.95 and 7 = 1.2. 



effective theory for the ripple modes is relatively valid for 
describing the imaginary part 5w(g) in Fig. 21 

The instability reflects the properties as a quantized 
vortex rather than those as a vortex sheet when the in- 
terface radius rg becomes small enough. In general, since 
the intercomponent interaction is essential to the shear- 
flow instability, it is expected that the instability is sup- 
pressed when the population of the second component 
is small (see i = 1 in Fig. [2|). Actually, the imaginary 
part of the numerical results is substantially smaller than 
that of the analytic estimation for L = 1 in Fig. SI The 
instability vanishes for i = 1 in the limit tq -^ for 
V — >■ i^'min without thc second component. Note that, for 
\L\ > 1, a dynamic instability can survive even in the 
limit To — >■ 0. The survived instability is well known as 
the vortex-splitting instability, which causes the splitting 
of a multiquantized vortex into single-quantized vortices 
[251 [26| . The abnormal behavior of the |/| = L mode 
above is regarded to be a reflection of the splitting insta- 
bility. In this sense, we call the |/| = L mode thc splitting 
mode in this paper. 

The instability of the core-flow vortex with L = 1 
is triggered by the vortex-characteristic modes: Kelvin 
(|/| = 1) and varicose {I = 0) modes. If the velocity of the 
second component or Q becomes smaller, only the Kelvin 
mode {\l\ = 1) has a nonzero imaginary part: 3w ^ 0. 
The effect of the Kelvin mode is qualitatively different 
from that of the ripple modes in the helical vortex sheet 
with larger i; although the ripple modes localized around 
r = rg do nothing other than cause ripples to propagate 
along the interface, the Kelvin mode makes the position 
of the core slide from r = 0. When Q is increased, the 
varicose mode (/ = 0) begins to have complex frequencies 
in addition to the Kelvin mode in the core-flow vortex 
with L = \. The varicose mode induces density per- 
turbations with rotational symmetry around the vortex 
axis in the density distribution | ^ j P . The perturbations 
may be considered as oscillations of the radius of the vor- 
tex core since the varicose mode causes the |^i| = |^2| 
plane to oscillate. We call this dynamic instability, trig- 
gered by the Kelvin and varicose modes in the core-flow 
vortex with i = 1, the core-flow instability (CFI) for the 
purpose of distinguishing it from the helical KHI in the 
helical vortex sheet regime. 

We summarize here the linear stability of the helical 
shear-flow states with a phase diagram in Fig. [51 The he- 
lical shear-flow states are dynamically stable below the 
critical value of Q or the critical velocity of the second 
component. It is surprising that the analytic result of 
the critical velocity (Vc) in Fig. |3]is in good agreement 
with the numerical result even in the crossover regime 
with small L > 1 in Fig. [51 This accidental agreement 
comes from the fact that the critical value (Q;=|l|) of the 
splitting mode is higher than that of the ripple modes, 
which is described well with the effective theory. Differ- 
ent modes have complex frequency for larger value of Q, 
leading to a more complicated instability. The simplest 
case is the core-flow vortex (L = 1), where the instability 
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FIG. 6: Phase diagram of the helical shear-flow instability 
for 7 = 1.2 and v — 0.95. The stable region is shown by 
filled circles . Different marks represents which modes can be 
amplified in the instability development. Filled diamonds are 
the analytical results in Eq. (|24p . 



is caused mainly by Kelvin or varicose modes. 



V. NONLINEAR DEVELOPMENT 

To gain deeper insights into the helical shear-flow in- 
stability, let us discuss the nonlinear development of the 
instability phenomena by numerically solving the GP 
equations (^5)1 . We restrict our investigation to the cases 
of the crossover or core- flow vortex regimes since the non- 
linear development in the helical KHI limit (r — >■ tq) must 
be essentially the same as that of the quantum KHI for 
a flat vortex sheet discussed in Refs. [3, H]- 

Thc numerical simulations were performed on a pe- 
riodic system along the z axis. Similar situations can 
be realized experimentally by considering, e.g., elongated 
condensates in a cigar-shaped potential. To neglect the 
flnite-size effect and capture the essence of the instabil- 
ity, we investigate the developments without a potential 
gradient around the z axis by using the cylinder potential 
V,- = Vtrap tanh[(r2 - Rl^p)/^^] with Ftrap - 50fii and 
-Rtrap ^ £,■ Unless otherwise noted, all the time devel- 
opment diagrams displayed in the following are obtained 
by using numerical simulations starting from stationary 



states $j with a small amount of white noise added to 
trigger the dynamic instability, which breaks the trans- 
lational and rotational symmetry along the z axis. The 
time evolution is obtained by solving the GP equation 
with the Grank-Nicolson method. The parameters of the 
initial states are flxed as 7 = 1.2 and v = 0.95. 

As a flrst step, we demonstrate the instability develop- 
ment from the core-flow vortex {L = 1). This instability 
corresponds to a GFI triggered by the Kelvin and varicose 
modes, which is distinguishable from the instabilities by 
the amplifications of thc splitting and the ripple modes. 
We discuss the instability for the case of L > 1 after the 
CFI. 



A. CFI by Kelvin mode amplification (L = 1) 
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FIG. 7: Nonlinear development of the Kelvin mode ampli- 
fication in the CFI of the core-flow vortex for L = 1 and 
Q^ = 0.98. The surface shows the isosurface |*ip = 1*2^ 
and the surface color represents the phase of the second com- 
ponent on the isosurface. The system size is 25.6^ in the x 
and y directions and the periodicity is 51.2^ in the z direction. 

The simplest case is the CFI developing from a core- 
flow vortex, where only the Kelvin modes have com- 
plex frequency. Figure [7] shows the typical develop- 
ment of the CFI owing to Kelvin mode amplification for 
Q^ = 0.98. The dynamic instability amplifies mainly the 
Kelvin mode with the largest value of the imaginary part 
5a; in this system, which is obtained by numerically solv- 
ing the BdG equation. The amplification leads to a de- 
formation of the initial straight vortex line [Fig.[7l^a)] into 
a helix [Fig. UV^)] with a periodicity that corresponds to 
the wavelength of the most amplified Kelvin mode. The 
radius of the helix grows with time and the growth stops 
up to a finite radius. After the growth, the core-fiow vor- 
tex typically keeps its helical structure although the helix 
is somewhat deformed [Fig.[71Jc)]. 

The direction of the helix appearing in the CFI is 
uniquely determined. To understand this fact qualita- 
tively, let us parameterize the trajectory s of the helix 
with z as 



s{z) = [rs cos{ksz),rsSm{ksz), z] 



(30) 



in Cartesian coordinates. Here, Ts and ks are the ra- 
dius and the wave number of the helix, respectively. It 
is important to note that the Kelvin mode amplification 
reduces the angular momentum of the first component 
containing a vortex because the amplification leads to a 
displacement of the vortex core from the z axis. During 
this process, the angular momentum of the second com- 
ponent must increase since the sum of the angular mo- 
menta of the two components is conserved in our system 
with rotational symmetry about the z axis. In the initial 
straight vortex state of Eq. (|4]) with Q > and L > 0, 
the first component has a positive angular momentum 
about the z axis. Since then the angular momentum of 
the second component becomes positive for kg > and 
negative for fc^ < 0, the helix must be right-handed with 
fcg > in order to conserve angular momentum. There- 
fore, whether the helix becomes right-handed (fcg > 0) or 
left-handed {kg < 0) is determined by the signs of Q and 



L; e.g., kg < for the core-flow vortex with Q < and 
L>0. 

Since there remains a relative velocity along the vor- 
tex line after the helical deformation, one may expect 
the CFI to occur additionally at the local point along 
the core-flow vortex. But this is not the case. Note that, 
in Fig. [71 the winding number of the phase arg ^2 of the 
second component displayed is kept through the insta- 
bility dynamics. In that case, the current velocity of the 
second component parallel to the helical trajectory is es- 
timated as ^Q/-\/rT"fc|rf. Thus, the relative velocity 
along the vortex core decreases with the radius r^ of the 
helix. Consequently, the helical deformation works to 
suppress the additional instability caused by the relative 
velocity along the vortex core. 

The growth dynamics of a helical vortex line in the 
CFI is similar to that in the Kelvin wave instability [ij] . 
However, the CFI is essentially a different phenomenon 
from the Kelvin wave instability. The CFI belongs to 
a dynamic instability caused by an internal interaction, 
the intercomponent interaction between the two compo- 
nents, and thus, the sums of the momenta and the angu- 
lar momenta along the z axis of the two are, respectively, 
conserved in our system with translational and rotational 
symmetry about the axis. In contrast, those quantities 
are not conserved in the Kelvin wave instability, where 
the thermodynamic instability is caused by energy dissi- 
pation resulting from external interaction with the envi- 
ronment. This essential difference between the two insta- 
bilities appears in their typical nonlinear developments; 
the radius of the helix stops growing within a finite range 
in the CFI, but, in the Kelvin wave instability, the ra- 
dius continues to increase monotonically unless the vor- 
tex collides with other vortices or the system reaches a 
local minimum of energy. 



B. Influence of varicose motion (L = 1) 

When Q is increased, varicose modes begin to have 
complex frequencies. Figure [5] shows the nonlinear devel- 
opment of the CFI for Q^ — 1.23, where both Kelvin and 
varicose modes have complex frequencies. In addition to 
the helical deformation of the vortex core resulting from 
Kelvin mode amplification [Fig.ISJb)], varicose mode am- 
plification causes oscillations of the radius of the vortex 
core [Fig. [Sl^c)]. This oscillation comes mainly from den- 
sity modulation of the second component trapped along 
the vortex core. We found that the winding number of 
the phase arg ^2 across the system in the z direction 
changes from Fig. [5{a) to Fig. IH^d). After the relative 
velocity decreases because of the reduction of the wind- 
ing number, oscillations of the core radius are essentially 
suppressed but the helical deformation of the vortex core 
survives. 

To better understand how the varicose mode amplifi- 
cation causes the reduction of the winding number in the 
second component, we examined the instability develop- 
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FIG. 8: Nonlinear development of the CFI induced by Kelvin 
and varicose modes for L = 1 and Q^ = 1.23. The system 
size is 25. 6^ in the x and y directions and the periodicity is 
51.2^ in the z direction. 



ment starting from the stationary state with an initial 
perturbation, which does not break rotational symme- 
try around the z axis. This perturbation results purely 
in the amplification of varicose modes without Kelvin 
mode amplification keeping the initial rotational symme- 
try of the condensate densities. Figure [9] demonstrates 
the nonlinear development of the second component in 
the CFI induced only by the varicose mode amplifica- 
tion for Q^ ~ \A1 . The density modulation caused by 
the varicose mode amplification makes deep troughs in 
the density of the second component. The phase arg^2 
drastically changes around the density troughs, and then, 
the winding number of the phase across the system along 
the core-flow vortex is decreased after an annihilation of 
a vortex-antivortex pair in the cross section of the phase 
distribution in Fig. [HI This phenomenon is an analog of 
the phase slippage of a superflow in a narrow channel, 
where the winding number of the superflow along the 
channel decays after vortices pass across the channel. 
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FIG. 9: Time development of the second component in the 
CFI from varicose modes for L = 1 and QE, = 1.47. Top and 
bottom panels show the close-up longitudinal profiles (y — 0) 
of the density |'I'2p and the phase 62 — arg»I/2 along the 
vortex core, (a) The varicose mode amplification leads to 
density troughs in the density profile along a narrow channel 
made by the vortex of the first component. The phase gra- 
dient is steep in the density troughs, (b) Two vortices from 
the outside approach each other around the trough as seen 
in the phase O2. (c) After the vortices annihilate, the phase 
becomes smooth and the relative velocity between the two 
components is decreased. 



C. Crossover regime (L > 1) 

For L > 1, the ampliflcation of splitting modes (|/| = 
L) with complex frequencies can have a dominant effect 



on the dynamic instability. In particular, as was men- 
tioned above, the splitting modes survive even when the 
second component is absent in the limit r ^ tq. Because 
we are interested in the phenomena characteristic not 
of the multiquantized vortices but of helical shear flows, 
we demonstrate here the instability dynamics where the 
splitting mode amplification is subdominant in the heli- 
cal vortex sheet regime. 
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FIG. 10: Development of the shear-flow instability for L — 2 
and Q^ = 0.74. Top panels show the density isosurface 
1^1 1 = |'I'2| and the phase Q2 = arg»I'2 of the second com- 
ponent on the isosurface. Bottom panels represent the distri- 
bution of the density difference ["I'll — |'I'2| around r = in 
the z — cross section. A circle in the cross section corre- 
sponds to the region containing vortices in the phase arg^i 
(not shown), (a) The cylindrical interface is straight in the 
initial state, (b) Amplification of the |/| — 1 mode causes the 
helical deformation of the interface, (c) Vortex splitting is vis- 
ible in the phase arg^i but not in the density difference, (d) 
The cylindrical interface partly splits into two singly quan- 
tized core-flow vortices in the place where the population of 
the second component becomes small locally, (e) Finally, the 
two core-flow vortices are twisted around each other. The 
system size is 40^ in the x and y directions and 51. 2^ in the 
z direction. 

Figure [TO] shows the dynamic instability developing 
from the helical shear flow for L = 2 and Q^ = 0.74. The 
ampliflcations of three types of modes (|/| = 0, 1,2) are 
coexistent in the dynamics according to the plot of L = 2 
in Fig. m A |Z| = 1 mode has the largest amplification 
rate in the state, and it induces a helical deformation of 
the cylindrical interface. In the early stage of the devel- 
opment, a deformation of the interface resulting from the 
ampliflcation of the |?| = 1 mode dominates over those of 
other modes [Figs. ITOTb) and 10(c)]. The radius of the 
cylindrical interface becomes thin locally as a result of 
density modulation caused by the varicose mode amplifi- 
cation. Then, vortex splitting occurs from the thin region 
[Fig.ITOId)]. In the final stage [Fig.ITOle)], a double helix 
of core-flow vortices is formed, and the varicose motion is 
suppressed by the reduction of the relative velocity along 
the helical trajectory of the core, as is the case with the 
instability of the L ~ 1 core-flow vortex in Fig. [5] 

The time development of a L = 3 core-flow vortex 
becomes more complex than that of the L = 2 core- 
flow vortex. Although the process of the development 
depends on which modes are dominantly amplified in the 
first stage of the instability, we essentially obtained the 



same result in the final stage: A triple helix of L = 1 
core-flow vortices was formed. 

The instability dynamics should change to those of 
quantum KHI when the interface radius tq increases with 
L. For our case, the signature of quantum KHI appears 
for L = 4. Figure [Tl] shows the instability development 
of the helical shear flow for L = 4 and Q^ = 0.74. In the 
early stage, a characteristic wavy pattern of quantum 
KHI [3, |j] appears at the interface in the cross section of 
the cylindrical interface in Fig. lTlT c). Since an |/| = 2 rip- 
ple mode has the largest amplification rate, there appear 
two waves in the interface pattern and an characteris- 
tic twisted structure is formed in three dimensions. Two 
single-quantized vortices are nucleated outside from the 
tops of the waves [bottom of Fig. [TIT d)]. and then, the 
nucleated vortices helically wind the two remaining vor- 
tices around the z axis [top of Fig. lTTT d')]. The remaining 
vortices form two single-quantized core-flow vortices and 
finally a helically twisted bundle of four core-flow vortices 
appears [Fig.fTlTe)]. The population of the second com- 
ponent is not distributed equally among the four vortices 
and the density modulation is so strong that the second 
component sometimes forms a large droplet in the cores 
along some vortices. 
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FIG. 11: Development of the instability from the helical vor- 
tex sheet for L = 4 and Q^ = 0.74. (a)-(c) The straight cylin- 
drical interface in the initial state is deformed into a unique 
three-dimensional structure as a result of the amplification 
of the 1^1=2 mode. In the 2 = cross section, the interface 
forms a wavy pattern characteristic in quantum KHI. (d) Two 
single-quantized vortices are released outside from the tops of 
the waves, and then the nucleated vortices helically wind the 
two vortices, (e) Finally, a twisted bundle of four core-fiow 
vortices appears with strong density modulation in the sec- 
ond component. The system is 76.8^ in the x and y directions 
and, 51. 2^ in the the z direction. 



VI. CONCLUSIONS AND DISCUSSION 

The helical shear flows can be classified into two states: 
a helical vortex sheet, where the radius of a cylindrical 
interface between the two components is large and the 
vorticity is localized at the interface, and a core-flow vor- 
tex, where the population of the inside component flow- 
ing along the core becomes smaller and the vorticity is 
distributed around the core. These states develop into 



exotic instability phenomena in which various kinds of 
modes compete. 

The linear stability of helical vortex sheets is well de- 
scribed with the effective theory of the helical KHI. When 
the supcrfluid velocity of the inside component exceeds 
a critical value, ripple modes trigger the dynamic insta- 
bility, called the helical KHI. When the interface radius 
becomes comparable with the thickness of the interface, 
our numerical analysis revealed that the instability phe- 
nomena are affected strongly by the vortex-characteristic 
modes: Kelvin, splitting, and varicose modes. 

We showed the nonlinear development of the shear-flow 
instability by using three-dimensional numerical simula- 
tions. If the outside component has a single-quantized 
circulation in the core-flow vortex regime, the instability 
transforms the core-flow vortex from the initial straight 
line into a helix due to amplification of Kelvin modes. 
The direction of the helix is uniquely determined de- 
pending on the flow directions of the inside and outside 
components. The helical deformation reduces the local 
relative velocity along the vortex core between the two 
components to suppress additional instabilities on the 
core-fiow vortex. When the velocity of the second compo- 
nent becomes larger, varicose modes in addition to Kelvin 
modes arc amplified in the instability development. The 
varicose modes induce density modulations, leading the 
phase slippage in the second component to decrease the 
relative velocity more effectively than the helical defor- 
mation. These instabilities of the core-fiow vortex caused 
by the amplification of Kelvin and/or varicose modes are 
called the CFI. In the crossover regime between the he- 
lical KHI and the CFI, the primary difference between 
the KHI and the CFI is caused by the splitting modes, in 
which multiquantized vortices split into single-quantized 
vortices. The vortex splitting effect leads to a helically 
twisted bundle of single-quantized core-fiow vortices in 
the final stage of the instability. 

The helical shear-fiow instability discussed in this pa- 
per is one of the fundamental hydrodynamic instabilities 
in multicomponent supcrfluid systems. This instability 
plays an important role in nonequilibrium dynamics of 
phase-separated two-component BECs containing quan- 
tized vortices when there is a larger population imbalance 
between the two components. For example, numerous 
core-fiow vortices are nucleated via an annihilation of two 
interfaces (a domain wall and an antidomain wall) (ij- 
\u\ . In particular, the CFI must be crucial to understand 
the dynamics and stability of a vorton [23, [2g , a loop of 
core-fiow vortex. These situations can be realized in two- 
component BECs with the experimental technique that 
was followed in Ref. [29[, where the nodal plane in one 
component was filled with the other component and then 
the filling component was selectively removed by using a 
resonant laser beam. Moreover, details of the instabil- 
ity developments can be investigated experimentally in 
cigar-shaped condensates, in which a multiquantized vor- 
tex with large circulations is prepared in a component by 
using Laguerre-Gaussian beams or topological phase im- 



10 



printing [30I433J . Relative motion along the vortex core 
can be realized, e.g., by employing a technique similar to 
that of Ref. Q by utilizing the different Zeeman shifts be- 
tween the two components. These realizations of the he- 
lical shear-flow instability provide an essential step prior 
to establishing modern quantum hydrodynamics 5| and 
lead to new physical insights not only into multicompo- 



nent BECs but also into other multicomponcnt supcrfluid 
systems. 
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